Entanglement Entropy of Systems with Spontaneously Broken 

Continuous Symmetry 



Max A. Metlitski^ and Tarun Grover^ 

^Kavli Institute for Theoretical Physics, 
University of California, Santa Barbara, CA 93106 
^Department of Physics, University of California, Berkeley, CA 94720 

(Dated: December 23, 2011) 



Abstract 

We study entanglement properties of systems with spontaneously broken continuous symmetry. 
We find that in addition to the expected area law behavior, the entanglement entropy contains a 
subleading contribution which diverges logarithmically with the system size in agreement with the 
Monte Carlo simulations of A. Kallin et. al. (arXiv:1107.2840). The coefficient of the logarithm is a 
universal number given simply by ^2~^^ where Nq is the number of Goldstone modes and d is the 
spatial dimension. This term is present even when the subsystem boundary is straight and contains 
no corners, and its origin lies in the interplay of Goldstone modes and restoration of symmetry in a 
finite volume. We also compute the "low-energy" part of the entanglement spectrum and show that 
it has the same characteristic "tower of states" form as the physical low-energy spectrum obtained 
when a system with spontaneously broken continuous symmetry is placed in a finite volume. 
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I. INTRODUCTION 



In recent years there has been a lot of theoretical interest in entanglement properties of 
quantum states of matter. Entanglement has proved to be a useful probe of non-local corre- 
lations for both gapless and gapped systems. The most commonly used characterization of 
entanglement is the entanglement entropy S — —ir(pA log p^), defined as the von Neumann 
entropy associated with the reduced density matrix pA of a subsystem. A close relative of 
the entanglement entropy is the Renyi entanglement entropy Sn — — log^^lp^)- 

To date, the most impressive progress has been made in the study of entanglement entropy 
of one dimensional critical states. Here, it has been shown^'^ that for systems described by 
conformal field theories (CFT's) the entanglement entropy of a system of length L behaves 
as. 



where a is the short distance cut-off. The coefficient c is a universal number known as the 
central charge of the CFT. The subleading constant 7 depends on the system geometry {e.g. 
the ratio of the subsystem size to the full system size). Although 7 is not fully universal 
as is clear from the cut-off dependence of (1.1), it is universal up to an additive constant. 
The universal behavior of the entanglement entropy (1.1) has proved useful for extraction of 
the central charge of the governing CFT in numerical density-matrix renormalization group 
(DMRG) studies of one-dimensional critical systems. 

Our present understanding of entanglement in dimension d > 1 is far less complete. 
However, it is generally expected that for both critical and non-critical systems the leading 
contribution to the entanglement entropy scales as the area of the subsystem boundary, 
S — CA/a'^~^.^ This "area law" contribution is related to short-range entanglement in the 
vicinity of the boundary, and as a result, the proportionality constant C is non-universal. 
However, for critical scale invariant systems one expects a subleading, fully universal, geo- 
metric contribution 7 to the entanglement entropy. 



The scaling form (1.2) is believed to hold for scale invariant systems in = 2 with arbi- 
trary smooth boundaries and for systems in d = 3 whose boundary is flat."''"^ Additional 
logarithmic contributions are expected in = 2 if the boundary has sharp corners and in 
d = 3if the boundary is curved. We note that all the above stated results/hypotheses on the 
entanglement entropy apply also to the Renyi entanglement entropy Sn, with the constants 
C and 7 acquiring a dependence on n. 

The scahng hypothesis (1.2) relies on the following argument. The entanglement entropy, 
being a dimensionless quantity, can only depend on ratios of length (or energy) scales. 



A notable exception is provided by systems with a Fermi-surface, where the area law receives a multi- 
plicative logarithmic correction. ^'^ 




(1.1) 



S = CA/a'^-^ + 7 



(1.2) 
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However, in a scale invariant theory, the only two length scales are the total system size L 
and the short distance cut-off a, with the corresponding energy scales and where 
z is the dynamical critical exponent. Thus, any dependence of the entanglement entropy on 
the system size must come together with the dependence on the short distance cut-off. The 
variation comes from short-range entanglement in the vicinity of the boundary and is 
expected to take the form of some local geometric quantity integrated over the boundary 
area. If the boundary is straight and has no corners, there are no non-trivial local geometric 
quantities, and the only possible cut-off dependent contribution to the entanglement entropy 
is proportional to the integral of 1 over the boundary, i.e. the boundary area, in accord 
with Eq. (1.2). On the other hand, for a two-dimensional boundary with corners each 
corner can contribute a constant to resulting in an entanglement entropy which depends 
logarithmically on a and hence on L. For a generalization of these arguments to curved 
boundaries, see Refs. 9-11. 

In a recent breakthrough it has become possible to extract the Renyi entanglement en- 
tropy of quantum systems using Monte-Carlo simulations.^^ One of the first applications of 
the method of Ref. 12 has been to study the entanglement entropy of a spin- 1/2 Heisenberg 
model on a two-dimensional square lattice. The Monte-Carlo simulations were performed 
using a torus geometry, with the subsystem being either a cylinder or a square. As expected, 
the leading contribution to the Renyi entanglement entropy was found to scale linearly with 
the system size. However, surprisingly, a subleading logarithmic correction was observed for 
both geometries studied. 

The ground state of the Heisenberg model on a two-dimensional square lattice is known 
to spontaneously break the SU{2) spin-rotation symmetry to a U{1) subgroup. Thus, in the 
infinite volume limit the ground state is infinitely degenerate and the ground state manifold 
is a two-dimensional sphere, whose points are labeled by the orientation of the Neel order 
parameter n. The low energy excitations above a particular ground state are two linearly 
dispersing Goldstone bosons, known as spin-waves. The interactions between the spin-waves 
are irrelevant in the low-energy limit and one, thus, simply obtains a theory of two free scalar 
bosons. Naively, this theory is scale invariant and so according to Eq. (1.2), one does not 
expect any subleading logarithmic contributions to the entanglement entropy as long as the 
boundary of the subregion is smooth, e.g. for the cylindrical subregion geometry. For the 
square subregion geometry, one could attribute the logarithmic correction to the corners of 
the boundary, however, the Monte-Carlo estimate of the coefficient of the logarithm differs 
both in sign and by one order of magnitude compared to the previous calculation based on 
a free bosonic field theory.^ 

The caveat to the above discussion is that a system with a spontaneously broken contin- 
uous symmetry is, in fact, not scale invariant at low energy in the same way that a CFT 
is. The reason for this is that symmetry is always restored in a finite volume and instead 
of the degenerate vacuum manifold one has a unique ground state and a "tower" of excited 
states carrying a definite charge under the symmetry group. For instance, in the case of the 
Heisenberg model with its SU{2) U{1) symmetry breaking, the tower of states can be 
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described by an effective Hamiltonian, 



(.2 5-2 

— * 

where 5" is the total spin of the system, c is the spin-wave velocity, ps is the spin-stiffness 
and V is the volume of the system. Note that the spacing between the energy levels of the 
tower scales as — L"*^. For system dimension d > 1, where spontaneous breaking of 
continuous symmetry is possible, this spacing is much smaller than the spin- wave gap c/L. 

Thus, in a system with spontaneous breaking of continuous symmetry, to a given length 
scale L there correspond two infra-red energy scales: the tower of states spacing c^/ {psL'^) 
and the spin- wave gap c/L. In general, we expect the universal part of the entanglement 
entropy AS to depend on the ratio of these energy scales: 

S^C^^+ASip,L'-'/c) (1.4) 

a 

Here, the function AS* a priori can depend on the system geometry. Note that the argument 
of AS* is much greater than one.^ In this paper, we argue that for a system with 0{N) 
0{N — 1) symmetry breaking, A^* behaves as 

A5 = Mog(p,L'^-Vc)+7ord, ^ = (1-5) 

As usual, Eq. (1.5) assumes smooth boundaries in d — 2 and flat boundaries in d — 3. Thus, 
the coefficient of the logarithmic divergence with system size L is directly expressed in 
terms of the number of Goldstone modes N — 1. In particular, for the case of the Heisenberg 
antiferromagnet N — 3, b — 1. We note that in addition to the logarithmic correction 
in Eq. (1.5) there also appears a finite constant 7ord that is generally expected to depend 
on the system geometry. Unlike in the case of entanglement entropy in one dimension, 
Eq. (1.1), where the presence of a logarithm rendered 7 universal only up to an additive 
contribution, here 7ord is fully universal, as all the short-distance physics is absorbed into 
the order-parameter stiffness Ps and the Goldstone velocity c. We will present an expression 
for 7ord in terms of a related constant in a free scalar theory. We numerically evaluate the 
latter constant (hence obtaining 7ord) using the correlation matrix technique of Ref. 14 for 
the geometry studied in the Monte-Carlo simulations of Ref. 13. 

In addition to the entanglement entropy, we compute the Renyi entropies Sn, which are 



^ Clearly, for L — >■ 00, the quantity psL'^~^/c ^ 1. One can, however, ask if this quantity can become of 
0(1) or smaller in a system close to a continuous phase transition into an antiferromagnetically disordered 
phase where ps — )■ 0. For the effective low-energy description of the antiferromagnet to apply, we need the 
system size L to be much larger than the correlation length of the system ^, thus, psL'^~^/c ^ ps^'^~^ /c. 
A priori, there need not be a relation between pg and ^. However, for all quantum phase transitions known 
to the authors, Ps^'^~^ /c either goes to a constant at the transition or diverges. 
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also found to satisfy the scaling form in Eqs. (1.4), (1.5). The coefficient of the area law C is 
now expected to depend on the replica index n and so is the finite size constant 7. However, 
the coefficient of the logarithmic correction h is found to be independent of the replica index. 

We also study the full low energy spectrum of the entanglement Hamiltonian He defined 
in terms of the reduced density matrix of the subsystem pA as pA = exp(— if^). We show 
that the low energy part of He up to a constant energy shift satisfies He = Htomer/TE-i with 
i^tower, Eq. (1.3), the physical tower of states Hamiltonian of the subsystem. The "effective 
temperature" Tg is found to scale with the system size as ~ c/L. These findings are 
remarkably similar to recent results on gapped topological phases, which possess gapless 
edge states. ^^"^^ It was found that for such phases the entanglement Hamiltonian associated 
with a subsystem is proportional to the actual Hamiltonian governing the edge states that 
would be exposed if a physical cut was performed along the subsystem boundary. 

Our conclusions are based on two calculations. The first involves a simple quantum 
mechanical model of two coupled quantum 0{N) rotors Ua and Ub, representing the average 
order parameter in each subsystem, 

„2 f 2 J2 f 2 

H-^ + ^-J^A-UB (1.6) 

Here La,b is the angular momentum of each rotor. Va, Vb and V denote the volumes of each 
subsystem and of the total system respectively. We choose the coupling J ~ PsL'^ '^ appro- 
priately to reflect the flnite order-parameter stiffness of the system. It is straightforward to 
flnd the spectrum of the Hamiltonian (1.6) in the limit psL^~^ jc ^ 1. The entanglement 
properties of the ground state wave-function are precisely as described above. 

Ideally, we would like to derive the effective Hamiltonian (1.6) starting from the 0{N) 
non-linear (j-model, which forms the full low-energy description of the system. At this 
point, we have not succeeded in doing so. Instead, we have performed a calculation of 
the entanglement entropy in the a-model itself using the replica method. For a general 
A^, taking the compact nature of the order parameter into account is non-trivial and our 
calculations are, strictly speaking, reliable only as long as the temperature T is much larger 

2 

than the tower of states energy spacing Ato^gr ~ -fy- In this limit, the replica method 
calculations of entanglement entropy are in agreement with the results of the rotor model 
(1.6). Moreover, the rotor model indicates that the leading behavior of entanglement entropy 
is temperature independent for T much smaller than the spin- wave gap c/L, so our replica 
calculations remain correct down to zero temperature. We confirm this fact explicitly for 
the case = 2, where replica calculations can be extended all the way to T = 0. For a 
general A^, we also study the entanglement entropy in the presence of a symmetry breaking 
field which lifts the degeneracy of the vacuum manifold, and find complete agreement 
between the rotor model and replica method results. 

We would hke to note that the presence of logarithmic corrections to the entanglement 
entropy of a Heisenberg antiferromagnet has been theoretically pointed out previously in 
Ref. 23. The authors of Ref. 23 have used the spin- wave (large S') expansion of the Heisen- 
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berg model. The effect of symmetry restoration in a finite volume has been mimicked by 
an application of a staggered magnetic field h ~ suitably chosen to give a zero net 

staggered moment. The final step of the calculation of the entanglement entropy has been 
performed numerically and finite size scaling was used to extract the coefficient of the loga- 
rithmic divergence b ^ 0.93. This value is quite close to our exact result b = 1.'^ It was also 
found numerically that the coefficient of the corresponding divergence in the Renyi entropy 
Sn is practically independent of n in agreement with our exact conclusion. 

The connection between the restoration of symmetry in a finite volume and the appear- 
ance of subleading logarithmic terms in the entanglement entropy has also been made in 
Ref. 13. Here the authors started with a mean-field Neel state of a Heisenberg antiferro- 
magnet and averaged it over all orientations of the order parameter. The entanglement 
entropy of the resulting spin-singlet state was found to be iS = log A/" = dlog{L/a), where 
A/" is the total number of lattice sites (here and below we drop the constant piece in S). 
In contrast, our calculation, Eq. (1.5), gives for a Heisenberg antiferromagnet {N = 3), 
AS ={d-l) log ((p,/c)i/("'-i)L). The physical reason for the difference is that Ref. 13 does 
not take into account the existence of spin-waves, while as we have argued above, the pres- 
ence of both the tower of states and spin- waves is crucial for the correction to entanglement 
entropy (1.4). The estimate of Ref. 13 for an antiferromagnet is physically similar to the 
exact result one obtains in the case of a free Bose gas where S — ^ log = | log L/a, with 

- the total number of particles and a - the average interparticle distance. On the other 
hand, for an interacting Bose g superfluid, our result, Eq. (1.5), with N = 2 gives 

S = log (^{ps/cY^^'^^^^ Ly Thus, based on entanglement entropy we can distinguish a free 
Bose gas, where no Goldstone boson is present, from an interacting superfiuid, which has a 
Goldstone mode. 

This paper is organized as follows. Section II presents a calculation of entanglement prop- 
erties of the toy rotor model (1.6). In section HI we perform a calculation of entanglement 
entropy in the 0{N) non- linear cr- model using the replica method and show that the result 
is in complete agreement with that of the rotor model (certain details of the calculation are 
given in appendix A). Some concluding remarks are presented in section IV. 

II. ROTOR MODEL. 

In this section, we calculate the Renyi entropy in the model (1.6), describing two quantum 
0{N) rotors and Ub of unit length. These rotors are taken to represent the average 
orientation of the order parameter in subsystems A and B. We choose the coupling J 



^ We have repeated the calculation in Ref. 23 for system sizes up to 100 x 200 sites and found that b 
approaches unity (within a percent) as the system size is increased. We also note that the calculation 
in Ref. 23 can be reformulated using the Schwinger boson representation of spin operators. This yields 

an SU{2) spin-symmetric mean- field ansatz wavcfunction for the antiferromagnetic state that breaks the 
SU (2) symmetry spontaneously in the thermodynamic limit. 
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between the rotors in the foUowing way. RecaU that when the system is placed in a box of 
size L, the energy cost to twist the order parameter by an angle 6 between the two sides of 
the box is related to the order-parameter stiffness Ps via, 



1 d'^E 



(2.1) 

61=0 



Thus, we take J ~ psL'^~'^. Note that this value of J is approximate and should be under- 
stood in a scaling sense only. 

To find the spectrum of the model (1.6) it is convenient to first work with the Lagrangian 
formulation, 

= 2^2 + 2^2 -^^A ■ ub (2.2) 

Let us introduce the average and relative coordinates n and 5na. via. 



riA 



(r) = n^/l - a^{6naY + aea{T)6na{T), a = Vb/ (Va + Vb) 



UBir) = n^/r^¥{5^ + bea{T)5naiT), b = -Va/{Va + Vb) (2.3) 

Here, — 1 and are {N — 1) unit vectors forming an orthonormal set together with 
n. We expect the fluctuations of the average coordinate n to be parametrically slower than 
those of the relative coordinate Sn. Moreover, we expect the fluctuations of 5n to be small. 
Therefore, expanding the Lagrangian (2.2) to leading order in 5n and in derivatives of n we 
obtain. 



Ln + Lsn (2.4) 

2c2 



Ln = ^{drnr (2.5) 



Lsn = ^{djn^y + "-{6n^y (2.6) 

Here V — Va + Vb is the total volume and = VaVb/{Va + Vb) is the reduced volume. We 
see that the Lagrangians of the average and relative coordinates decouple. The dynamics of 
the average coordinate n are governed by the quantum rotor Lagrangian (2.5); the associated 
Hamiltonian 

Hn — 2^ 'y ' -^totai — La + Lb, (2.7) 

is, indeed, the appropriate "tower of states" Hamiltonian describing the lowest energy exci- 
tations of a system with spontaneous symmetry breaking. On the other hand, the dynamics 
of the relative coordinate Sn are governed by the Lagrangian (2.6) describing an — 1 
dimensional harmonic oscillator. The frequency of this harmonic oscillator is given by 

J c 
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In the true physical system, the spectrum of the "relative motion" involves N — 1 Goldstone 

— * — * 

modes with a dispersion u = c\k\, where the momentum k is quantized in a finite geometry. 
Our rotor model (1.6) replaces this multi-mode spectrum with a single — 1 dimensional 
oscillator whose energy (2.8) is of order of the finite-size gap of the Goldstone modes. It 
turns out that such a replacement is sufficient for capturing the logarithmic contribution to 
the entanglement entropy (1.5). 

The ground-state wave-function corresponding to the Lagrangian (2.4) is a product of a 

— * 

-^totai — wave-function of the rotor Hamiltonian (2.7) and the ground state wave-function 
of the harmonic oscillator (2.6), 

^(^M - ^i=^(^^^5)^exp(-(<5nJV(2e^)) 

^ ^^^ (^^2)(l-i)/4 (-(^^ - ^Br/{2e)) (2.9) 
Here, 1-5^-^1 = 27r^/Vr(A^/2) is the volume of a unit sphere S^'^ and 

^2 \ 1/4 / ^ \ 1/2 



Note that the condition ^ <^ 1 guarantees that the amplitude of the relative fluctuations is, 
indeed, small. 

We proceed to compute the reduced density matrix from the wave- function (2.9), 



PA{nA,n'A) = J dnBi>{nA,fiB)'ip*{n'A,riB) (2.11) 



The contributions to the integral over Ub above come from — "^bI, \n'A — "^bI ~ ^ 1- 
Therefore, pA{nA,n'y^) is non-negligible only for ^ 1- We may change variables in 

Eq. (2.11) to 6n = Ha — ub and, to leading order, take both 6n and fiA — Ua to lie in the 
tangent plane of Ha- The integral over Sn then becomes Gaussian and gives, 

PA{fiA: u'a) = exp {-{Ha - fi'^)/{Ae)) (2-12) 

Since = — is just the Laplacian on the sphere, we may use the heat kernel expansion 

(n|e-^>0 ^ (4^,)iv-i)/2 exp(-(n-n07(4^)), s^O (2.13) 



Thus, as ^ <^ 1, we may write. 



(4^^2)(Ar-l)/2 

PA ^ ^ ^"^ " (2-14) 
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Defining the entanglement Hamiltonian He as pA — exp(— if^), we obtain, 

He = C^L^ + const (2.15) 

Hence, the entanglement Hamiltonian has the same form as the physical tower of states 
Hamiltonian. We may write He ~ Ha/Te where Ha — 2p Va Hamiltonian of the 

subsystem A and Tg — y^^2 ~ f ■ Therefore, the reduced density matrix (2.14) describes 
a thermal state of the subsystem with the temperature of order of the Goldstone gap c/L. 
We can now compute the Renyi entropy of the rotor model. Prom Eq. (2.14) 

C4 e2Nn(7V-l)/2 /4 f2N(n-l)(iV-l)/2 

triPA) - ^ 7^Li|n tr{e--' '^) ^ (2.16) 
where we've used Eq. (2.13) in the last step. Therefore, 

S„ = -;^logtr«) = ^ (log + ^) +log|5«-| (2.17) 

z c 

Note that if we wish to interpret the result (2.18) in terms of the actual physical system rather 
than the toy rotor model, we must remember that Eq. (2.18) only captures the contribution 
of the universal low-lying part of the entanglement spectrum given by Eq. (2.15). There will 
also be a contribution from the non-universal, high-lying part of the entanglement spectrum, 
corresponding to short-range entanglement in the system. This contribution is expected to 
give a standard area law term in the Renyi entropy, so we identify the Renyi entropy of the 
rotor model, Eq. (2.18), with the AS contribution in Eq. (1.4). In addition, since the result 
(2.17) depends on ^ (J), and hence on the details of the spectrum of "relative fluctuations" 
(Goldstone modes), it is only logarithmically accurate. In particular, it does not capture 
the universal geometric constant 7ord of Eq. (1.4). 

It is easy to generahze the calculation above to the case when the system is at a finite 
temperature T <^ ou ^ c/ L. In this case, the Hamiltonian of the relative coordinate Sn is 
still in its ground state, while the average coordinate ri is in a thermal state. Thus, the total 
density matrix of the system is, 

p{n,Sn;n',Sn') = l-{n\e-i^\n') x exp (-(((5nJ^ + (K)')/(2a) 



(2.19) 



where = tr(exp(— yr )) ^^'^ ^- ^' '^^^ obtained from nA,nB,n'j^,n'g by inverting 
Eq. (2.3). To calculate the reduced density matrix we must set hb = n's- The second factor 
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in Eq. (2.19) then implies that the density matrix vanishes for — nsl, \n'^ — "/^bI 

the other hand, the first factor in Eq. (2.19) varies on scale \n — n'\ 0(1) for T < c^/^psV) 

and on scale \n — n'\ ~ 

(cV(p,\/T))i/2 for T > c'^/ipsV). In either case this scale is much 
larger than ^ provided that T <^ c/L. Hence, we may set fiA = n'A = in the first factor, 

1 __£££L 1 

- \sL\ i^^2yN-m i-ii^A - nsf + {fi\ - fisf)/{2e)) (2.20) 

Note that the temperature dependence has disappeared in Eq. (2.20). Hence, for T <^ c/L 
the Renyi entropy is to leading order given by our T = result (2.18). One may understand 
this fact in the following way: by tracing over the subsystem B one puts the subsystem A at 
an effective temperature Te ~ c/L, therefore, an initial physical temperature of the overall 
system T <^Te can be neglected. 

Finally, let us consider the situation when the 0{N) symmetry of the system is explicitly 
broken by an external field h, which couples to the order parameters in the two subsystems 
as 

SL = -M{VAh-nA + VBh-nB) (2.21) 

where M is the expectation value of the microscopic order parameter {e.g. the staggered 
magnetization in the case when the system is an antiferromagnet). As before, changing 
variables to relative and average coordinates (2.3), we obtain, 

5L fa -MVh ■ n (2.22) 

Here, we assume that MVh <^ J such that to leading order the external field does not 
modify the dynamics of the relative coordinate. At the same time, we will take MVh 
to be much larger than the tower of states spacing c^/ (psV) so that the fiuctuations of 
the average coordinate n about the direction of the external field are small. Then writing 
fi — (tt, Vl — TT^) and expanding Eqs. (2.5), (2.22) in n we obtain, 

L ~ Ln + 

+ (2-23) 

with Lsn still given by Eq. (2.6). The Lagrangian (2.23) describes a,n N — 1 dimensional 
harmonic oscillator with frequency m given by 



Mc'^h 



1/2 

m = ( ^^^^ ) (2.24) 



Note that the condition c^/^psV) -C MVh <C J ensures that m is much smaller than 
the frequency of relative motion cu (2.8), but much larger than the tower of states spacing 
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/{psV). Now, the ground state wave- function of the system is 



(^^^^/(iv-i)/2 (-(1 - ^ • ^)/^^) e^P (-(^^ - ^B)V(2e')) (2.25) 



where 

.2 X 1/4 



c 



j (2.26) 



The conditions above ensure that ^ <^ <^ 1. Computing the reduced density matrix, 
Eq. (2.11), we observe that as ^ <^ we may replace n in the first exponential in Eq. (2.25), 
by ua- Then integrating over n^, 

pA{nA, u'a) ~ ^^^T^(N-i)/2 ®^P (-(2 - /i • (n^ + exp (-(riA - ^a)V(4C^)) 

^ (^^2)(l-i)/2 ^^P (-(^A + ^A)/(2eD) exp (-(Tf^ - ^'Afim) (2.27) 

where we've written ua = {t^Ai ^/l — 7?^) in the last step. Now, consider a harmonic oscil- 
lator Hamiltonian _^ 

He = ^ + \mEUjl^l (2.28) 
with the momentum conjugate to tta- For <^ 1, 

(Af-l)/2 

(7fA|exp(-i/^;)|7f^) j exp(-mEa;|(7fA + 7f^)V8)exp(-m£;(7fA-7rA)V2) 

~ (^)^ exp(-mEa;|(7fi + 7f^)/4)exp(-mE(7rA-7f;i)V2) 

(2.29) 

PA ~ a;;^~'e-^^ (2.30) 

12^ / X 

"^E = 777^, = (2.31) 

Note that ue <^ 1, permitting the expansion (2.29). Thus, we see that in the present case 
the low-lying part of the entanglement spectrum is that of a harmonic oscillator. However, 
unlike in the absence of an external field, the reduced density matrix is not a thermal density 
matrix of the subsystem. Indeed, the amplitude of oscillations of vf^ in the ground state 
of the entanglement Hamiltonian, (2.28), is = 1 / y/rriEUJE = V^^h, which is paramet- 
rically smaller than the physical amplitude of oscillations ^h- Hence, we cannot write the 
entanglement Hamiltonian in the form He = Ha/Te- 



Therefore, we may write 
with the identification. 
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It is now straight- forward to compute the Renyi entropy in the presence of an external 
field. From Eq. (2.30), 




(2.32) 



Hence, 

= log trp\ = (AT - 1) flog I + 

iV- 1 , c 

~ log 

2 ^mL 

where we've expressed the final result in terms of the physical energy gap of the total 
system m. Note that Eq. (2.34) crosses over to our zero-field result (2.18) when m becomes 
of the order of the tower of states spacing m ~ -fjji- As before, Eq. (2.34) should be 
interpreted as the contribution AS* in Eq. (1.4) to the Renyi entropy from the low-lying part 
of the entanglement spectrum. In addition, Eq. (2.33) is again by itself only logarithmically 
accurate since it involves the quantity ^. However, we observe that ^ cancels out in the 
difference of zero and finite field Renyi entropies, Eqs. (2.17), (2.33), 

s^^-si . (,„,_L _ i-l).io,|.-'| ^ ^ (lo,^ - 

(2.35) 

As we see, the difference of Renyi entropies (2.35) depends only on the physical quantities: 
the system volume, the order-parameter stiffness and the finite-field mass gap and, hence, 
is fully universal, capturing not only the logarithmic dependence but also finite terms accu- 
rately. In section III we will reproduce Eq. (2.35) by a calculation in the 0{N) non-linear 
(7-model. This agreement makes us confident that the simple rotor model presented in this 
section accurately captures the universal entanglement properties of the system. Moreover, 
we will interpret the Renyi entropy of a system in a field as that of a free-bosonic scalar field 
theory, providing a way to calculate the universal constant 7ord in Eq. (1.5) from a related 
constant in a free theory. 



(2.33) 
(2.34) 



III. REPLICA METHOD. 

In this section, wc calculate the entanglement entropy in the 0{N) non-linear cr- model 
using the replica method. The action of the theory is given by, 

S = ^J (fxdT {^{drfif + {ynf^ (3.1) 

with n- an A'"-dimensional unit vector. We will set the spin- wave velocity c — 1 below and 
restore it in the final results. 
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We choose the system geometry to be a rf-dimcnsional torus T'^. For simphcity, we take 
aU the directions of the torus to have length L. Wc choose the subsystem A to be the 
cyhndrical region i < x < L, with x - one of the directions of the torus. This is one of the 
geometries considered in the Monte-Carlo simulations of Ref. 13. Unless otherwise noted, 
we assume the ratio £/L to be finite. 

As is well known, the Renyi entropy Sn is given by 

Sn^ ^ log trp^^ = ^ log ^ (3.2) 

n — 1 n — 1 

where Zn is the partition function of the system on an n-sheeted Riemann surface. ^'^ For 
the geometry considered here this surface is given by (r, x) e (0, n^) x with /3 — 1/T - 
the inverse temperature. The following identifications need to be made on this space: 

{k(3^,x) {{k + l)l3',x), 0<x<e 
(0+,f) - (n/3-,f), £<x<L (3.3) 

Above, 0<A;<n — lisan integer. Below, we will be interested in the low temperature 
fimit ^ > L. 

To analyze the theory (3.1), we expand n about a given point on the unit sphere. 



n 



= (tt/VpI, V1-^Vp.) (3-4) 



with TT - an iV — 1 component real scalar field. At lowest order in energy expansion, the 
action (3.1) then becomes a free theory of the Goldstone modes tt, 

S^^J d'^xdT{d^Tif (3.5) 

Here, ^ runs over both space and time indices. In writing Eq. (3.5), we have lost the 
information about the compact nature of the order parameter manifold. We will partially 
restore the compactness later in the calculation. Thus, 

TV — 1 

log Z^ = —tr \og{-dX (3.6) 

where ((?^)n is the Laplacian on the n-sheeted Riemann surface. The calculation can be 
simplified by utilizing the translational invariance along the subsystem boundary, 

log Z„ = log((-^i)n + ^) (3-7) 

with k\\ = ^m, m G Z"^"^ the momentum along the boundary. The operator {d']^n now acts 
only in the two-dimensional (r, x) space. 
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So far, our calculation is identical to that in a free massless bosonic field theory. As 
long as the boundary conditions on the field are such that the energy gap is of order 1/L 
then one expects an entanglement entropy of the form (1.2). This was explicitly checked 
in Ref. 10, where twisted boundary conditions were assumed along the directions parallel 
to the subregion boundary, so that |A:||| > |;, with 6 -the twist angle. It is, therefore, 
clear that any violation of the scaling form (1.2) has to come from "quasi- zero" modes with 
eigenvalue of (— c?^)n much smaller than 1/L^. This is in accord with our discussion in 
the introduction, which concluded that the violation of the scahng form (1.2) in a state 
with a spontaneously broken continuous symmetry is due to an emergence of the tower of 
states energy scale Atower ~ which is much smaller than the spin- wave gap j^. Thus, 
we concentrate below on the contribution of quasi-zero modes to the entanglement entropy. 
As already noted, in a free theory quasi-zero modes can be eliminated by suitable boundary 
conditions, e.g. twisted boundary conditions on a torus or Dirichlet boundary conditions 
for open boundaries. However, for a physical system, in the absence of any external fields, 
the boundary conditions on n and, therefore, tt, will be periodic (non-twisted) along any 
periodic direction and von-Neumann for any open boundaries. Such boundary conditions 
admit a true zero mode with a constant tt , which will be discussed in more detail below. As 
we will now see, small temporal fluctuations on top of this zero mode give rise to quasi-zero 
modes. 

— * 

Let us compute the spectrum of quasi-zero modes. Clearly, these must have A;|| = 0. 
Let us now focus on the behavior of quasi- zero modes in the (r, x) plane. Here, we must 
deal with the branch cuts at r = A;/?, < x < The behavior of the eigenmodes in 
the vicinity of the branch cuts is expected to be nontrivial. However, for \t — k(5\ ^ L, 
we expect the eigenfunctions of — (9^)„ to approach a linear superposition of plane wave 
states, (f) ~ Q^'^r+2mmyij L ^ m G Z. For the eigenvalue — uP' -\- {^-nmjVf to be much smaller 
than l/L^, we must choose m = 0. Therefore, we expect the following asymptotic behavior 
of quasi-zero modes with eigenvalue — i 

0(t, x) = A+e''"(^-'='3) + A-e-*'"("-'^) , r - > L, (A; + l)/3 - r > L (3.8) 

Without loss of generality, we take oj > 0. To study the low temperature hmit /? S> L it is 
convenient to cut our Riemann surface into n separate sheets by defining 

0fc(T,x) = 0(A;^ + r,x) (3.9) 

Here the variable k is defined modulo n. Each sheet has a branch-cut atr = 0, 0<x<£ 
and the sheets are glued together along these branch cuts, 

0fe(O+,x) = (/>fc+i(0-,x), 0<x<£ 

a^(/.fc(o+,x) = a^(/.fc+i(o-,x), < x < £ (3.10) 
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At finite ^, we should also glue the sheets via 



4>k{T,X) = 4>k+l{T - /3,x) (3.11) 

However, since for |r| ^ L, approaches the plane wave states (3.8), we will take the r 
coordinate in each sheet to run from — oo to oo and implement Eq. (3.11) via a boundary 
condition, 

0,,(r, x) = A+e*'^^ + Ale-'^\ r ^ oo 

0fe(T, x) = AUe'-^^^^^ + A-_,e-^-(-+«, r ^ -oo (3.12) 

The corrections to this approximation are expected to be of order e^^^. Thus, we must solve 

= (3.13) 

subject to the boundary conditions (3.10), (3.12). 

We may further simplify Eq. (3.13) in the quasi-zero mode limit w <^ 1/L. Indeed, in the 
vicinity of the branch points, we expect to vary on the length-scale L. Thus, the typical 
contributions to the left hand side of Eq. (3.13) are of order 1/LF' and we may set the right 
hand side of Eq. (3.13) to zero, 

-9'0fc = O (3.14) 

We expect the approximation (3.14) to hold as long a |r| <C l/cj. On the other hand, once 
|r| 3> L the asymptotic forms (3.12) become vahd, so we should solve Eq. (3.14) subject to 
the boundary conditions, 

(t>k{T, x) = {At + A-) + i{Al - A-;)ujT, r ^ oo 

0,(t,x) = At,e-^ + A-_,e--^ + i(A^,e-''-A-_,e-^^^^ r ^ -oo (3.15) 

together with the gluing condition (3.10). 

We can solve the Laplace equation (3.14) using conformal mapping. Let us define z = 
x+ir on each sheet of the Ricmann surface. Topologically, each sheet is a cylinder due to the 
periodicity of x. We begin by mapping each cylinder to a complex plane using w — e^'^'^^l^ , 
This maps the branch cut atr^O, 0<x<£toa segment of a unit circle with angle 
< ^ < e^'^'^/^. Moreover, r = oo is mapped to the origin and r = — oo to the point at 
infinity. Next, we apply the following map to each sheet, 

\ — w 

This maps the branch cut into the positive x-axis. Moreover, r = ±oo is now mapped to 
_g±7ri^/L^ Finally, we map the n-sheets into a single sheet using s = C^^"- More precisely, 
for C on the A;-th sheet with ( = re'^, r > 0, < 6^ < 27r, we define C^^" = r^/"e^(^+2'''=)/". 
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With this definition and the gluing (3.10), there are no branch cuts in the s plane. Note 
that the points r = ±00 on the kth sheet now map into = e'^«(2fc+i±^/i)/n-^ 'jj-^g boundary 
conditions (3.15) imply logarithmic divergences for s — )■ s^. Note that these are the only 
singularities that can occur in the s-plane, including the point s = 00. Indeed, s = 00 is the 
image of the branch point r = 0, x = at which we expect no singularity. Therefore, 

= E (^fe 1^ - + ^fc log \'-'k\)+ C, (3.17) 

k 

Note that the absence of a singularity at s = 00 implies, 

E(^^ + ^fc")=0 (3.18) 

k 

The solution (3.17) has the following asymptotic behavior as r — > ±00 on kth sheet, 

0fc(T, x) ™ ~ + ( ^ sin(7r£/L) j 4 + ^^4 log I Sfe - I + E log | s+ - s^" | + C, 

Mr, x) + log (^sm{7ri/LU + log - sfl + E l^fe - «rl + ^ 

(3.19) 

Imposing the boundary conditions (3.15) we obtain, 

4 - ^(At-A-) (3.20) 



c 



k 



iuL 



(.Iti^^'^^ - A'k-ie-'-^) (3-21) 



and 



^fe -C = log(^sin(7r£/L))c+ + Eci+log|4-^/^l+E^rlog|4-^rl 

^ i^k I 

AUe'-^ + ^,-_ie— ^ - C = log (- sin(7r£/L)') c," + E 4 log - s+| + E log K - si\ 

^ I i^k 

(3.22) 

Note that from (3.21) the right-hand side of Eq. (3.22) is suppressed compared to the left- 
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hand side by a factor of a;L, and at leading order may be dropped. Thus, we must solve 



At + A- = C (3.23) 
A+_^e''"^ + ^^1^"'"^ = C (3.24) 
(1 - E - (1 - ^"'"') E ^fc" = (3.25) 

k k 

Here, Eq. (3.25) comes from the condition (3.18). To further simplify the above equations, 
we observe that our original problem is symmetric under the time translation r — >■ r + /9, 
which generates a cyclic permutation of the n sheets. In the s plane it maps s — )■ e^'^*/"^. 
We choose to be an eigenstate of this symmetry with eigenvalue e^'^*^'/", which implies. 

At = ^±6^'^^^'=/" (3.26) 

For p ^ 0, the symmetry also constrains C = 0. Moreover, in this case Eq. (3.25) is 
trivially satisfied. Then, solving Eqs. (3. 23), (3. 24) we obtain A+ — —A^ and the condition 
sin(a;/3) — 0, i.e. 

0; = ^, meN, p^O (3.27) 

Thus, for "momentum" p ^ 0, one mode is present for each m in Eq. (3.27). On the other 
hand, for p = 0, we obtain C = A_^_ + A_ and e^^ — 1, i.e., 

a; = ^^,meN, p = (3.28) 

In this case, A^ and A_ are independent, so two modes are present for each m in Eq. (3.28). 
The explicit form of the corresponding eigenstates is. 



A^ + A--'-^{A^-A-)J2^og 



27r 

k 



^-4 



^ A+ + A- + iujr{A+ - A-) (3.29) 



We have inverted the conformal mapping in the last step. Upon continuing Eq. (3.29) from 
\ujt\ <^ 1 to the full range of r we obtain, 

(f) = A+e'""^ + A-e-'"'^ (3.30) 

For u! satisfying (3.28), Eq. (3.30) is actually the exact eigenstate of —9^ on the n-sheeted 
surface. In contrast, the solutions with p 7^ that we've constructed are not exact. The 

leading correction to these eigenstates can be calculated by keeping the terms on the right- 
hand side of Eqs. (3.22). This gives an 0{L/P) corrections to the eigenvalues (3.27). We 
have verified that these corrections do not modify any of the leading order results below for 
the Renyi entropy. 

We are now ready to compute the contribution of zero and quasi-zero modes to the 
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partition function. 

logZ„ = log^o - ^^tr' \og{-d\ (3.31) 

As already mentioned, in addition to the quasi-zero modes discussed above, —d^ has one 
true zero mode on the n-sheeted Riemann surface, which is simply given by constant (p. Sq 
in Eq. (3.31) is the zero-mode contribution and tr' is the trace with the zero-mode removed. 
The trace in Eq. (3.31) requires regularization. We will use Pauh-Villars regularization 
here by subtracting the free energy on an n-sheeted Riemann surface of a free boson field 
with mass A, which acts as a cut-off. Since we are considering only modes with eigenvalue 
< l/L^, the cut-off A should be taken to be A ~ 1/L. Thus, 

logZ„ = log^o - ^^tr' log{-d\ + ^:i^tr log{{-d\ + A') (3.32) 

Note that the regulator trace includes the zero-mode of —9^. Let us begin by calculating 
the zero mode contribution Sq. Physically, the presence of the zero-mode indicates the 
degeneracy of the vacuum manifold. Indeed, we have expanded our field n around a fixed 
direction. However, in reality, we should integrate over all directions of n. Thus, 5*0 should 
be converted into an integral over the order parameter manifold. Locally the manifold can 
be parameterized by 

7r(x) = TTo (3.33) 



and 



with the metric 



9ij 



So^ J d^-^TToy/d^ (3.34) 



Thus, 



(3.36) 

where we have gone from a local integral over ttq to a global integral over the order parameter 
orientation n. Here, jS*^"^! = 271^^"^ /T(N/2) is the volume of an A?^ — 1 dimensional unit 
sphere. 
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Next, consider the contribution of quasi-zero modes, 

A^- 1 



2 

A^- 1 



{tr'\og{-dl)-tr'\og{{-d% + A')) 

oo 

- 1) 5^ (log(7rm//3)' - log + A')) 



[n 



m=l 



+ 2 J2 (log(27rm//3)2 - log {{27rm/pf + A^)) 



m=l 



(3.37) 



The first sum in the square brackets gives the contribution of modes (3.27) with p ^ 0, 
while the second sum gives the contribution of modes (3.28) with p — 0. Using the standard 
relation, 

^ oo 

-- {log{{27rm/f3f + uj')- log{{27rm/ (3 f + cu'^)^ log Zho{uJ,f3)- log Zhoico', (3) 

m=— oo 

(3.38) 

where Zhoi^^, is the partition function of a harmonic oscillator with frequency ou at inverse 
temperature (3, 

(3.39) 



1 



'ho 



2sinh(/3cj/2) 



we obtain in the limit /3A 1, 

N -1 



log^r 



-n/3A + (n + 1) log(/5A) + (n - 1) log 2) 



(3.40) 



Combining (3.40) with (3.36) and taking the contribution of the zero mode to the regulator 
trace in Eq. (3.32) into account, we obtain 



log^n 



N-1 



[n/3A + log(np,LV(27r^)) -(n-1) log(2^A)] + log (3.41) 



which gives the Renyi entropy. 



A5„ 



1 , Z„ 
log 



N-1 



n-1 ° Z" 
— — log 



log(p,L'^A/7r) - 



logn 
n-1 



+ log\S 



N-li 



(3.42) 
(3.43) 



Note that we have set the cut-off A ~ 1/L in the last step, so that Eq. (3.43) is valid only to 
logarithmic accuracy and we have dropped the constant terms in Eq. (3.42). We will come 
back to give an interpretation of these terms below. We have also restored the spin-wave 
velocity c in Eq. (3.43). 

A few comments are in order here. First of all, the result (3.43) takes into account only 
the contribution of zero and quasi-zero modes to the entanglement entropy. As argued in the 
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introduction, the higher energy modes are expected to give a CFT hke contribution (1.2). 
Second, the result (3.43) agrees with the expression (2.18) obtained in the toy rotor model 
of section 11. Third, our calculation above is strictly justified only for temperature much 
greater than the tower of states energy spacing, T ^ Atowcr ~ This can be seen from 

studying the physical free energy F{T) = — 4logZi. Observe that Eq. (3.41) gives, 



This is in agreement with the free energy of the rotor Hamiltonian (1.3) for T ^ Atowcr, but 
not for T <C Atower, where the physical F{T) — F{0) is exponentially suppressed. The tech- 
nical reason for the disagreement between our path-integral treatment and the Hamiltonian 
calculation is that the former takes the compactness of the order parameter into account 
only in the treatment of zero modes, whereas the finite frequency modes are still described 
by a non-compact free boson theory. This approximation is reliable for T ^ Atower where 
the fluctuations of n corresponding to the finite frequency modes are small, but breaks down 
for T below Atower where these fluctuations arc large. However, we observe that while our 
free energy has an unphysical temperature dependence for T ^ Atowcr > this temperature 
dependence disappears in the expression for the Renyi entropy (3.43). Thus, one may hope 
that the result (3.43) is correct not only for Atowcr ^ T ^ c/L, but actually for all T ^ c/L. 
The calculations of section 11 in the toy rotor model support this hypothesis. In appendix 
A we confirm this hypothesis by a replica method calculation for the special case N — 2, 
where the compactness of the order parameter manifold can be taken into account exactly 
for all temperatures. 

Let us also analyze the entanglement entropy in the presence of a small symmetry breaking 

— * 

field h coupled to the order parameter as. 



Here, M is the expectation value of the order parameter. In the case of a Heisenberg 
antiferromagnet the symmetry breaking field corresponds to a staggered magnetic field. 
Using the representation (3.4) and expanding the action to leading order in tt 



with the mass = hM/pg. As in section H, we will consider the regime when the field 
induced Goldstone mass is much smaller than the spin- wave gap 1/L, but much larger than 
the tower of states gap Atower ~ ^/{PsL'^)- In this regime, the spectrum of spin- waves with 
finite momentum is unaffected by the symmetry breaking field. On the other hand, the 
zero-momentum fiuctuations of the order parameter about the staggered field are now small 
and described by a harmonic oscillator with frequency m. Thus, in this field range, one 




(3.44) 




(3.45) 




(3.46) 
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may ignore the compactness of the order parameter manifold and work with the free theory 
(3.46). The regularized partition function on the n-sheeted Riemann surface now becomes 



logZ„ = {tr log{{-d\ + m') - tr log{{-d\ + A')) (3.47) 

We again consider only the contribution of zero and quasi- zero modes of (— c?^)„ to the 
partition function (3.47). Using the spectrum (3.27), (3.28) 



, r. N -1 

logZ„ = - 



(log((7rm//3)' + m^) - log ((7rm//3)' + A')) 



[n 

m=l 



2 

oo 

+ (log((27rm//3)2 + m^) - log {{27rm//3f + A')) + log(mVA') 



m=l 



(3.48) 



where the last term in the second hue of Eq. (3.48) denotes the contribution of the zero 
mode of (— 9^)„. Performing the sum and taking the limit T — > 0, 

- 1 

log = — — {n/3{A - m) - (n - 1) log A/m) (3.49) 
which gives the Renyi entropy, 

N -\ 

ASn = log A/m (3.50) 

Here, we have set A ~ 1/L in Eq. (3.51). Observe that Eq. (3.52) is in agreement with the 
result (2.34) of the rotor model. Also note that Eq. (3.51) implies that the Renyi entropy of a 
free theory diverges logarithmically in the massless limit m — > 0. In the case of a system with 
a spontaneously broken symmetry, the calculations above are, however, valid only as long as 
m is much larger than the tower of states gap Atower ~ j\d, and at m ~ Atower, Eq. (3.52) 
crosses over to our previous result (3.43). Since the present calculation is performed at 
T = it supports the hypothesis that Eq. (3.43) remains correct down to T = 0. 

Let us discuss the connection between the entanglement entropy of a free theory and a 
state with a spontaneously broken continuous symmetry in more detail. For a free single 
component real scalar field with mass m, we expect the Renyi entropy to behave as, 

Sn = + A,5free(mL) (3.53) 
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FIG. 1: The universal constant contribution 7free) Eq. (3.54), to the Renyi entropy 52 of a free 
massive scalar in two dimensions as a function of the ratio i/L. The total system is a torus with 
dimensions L x L, while the subsystem is a cylinder of size i x L. We take the mass m = 10~^ 
and study systems with size L up to 100 sites, so as to fit the universal part of the entanglement 
entropy to the scaling form in Eq. (3.54). 7f1.ee is related to the universal constant 7ord in the 
ordered phase via Eq. (3.55). 

with ASfree & Universal geometry dependent function.^ Here we assume all the usual condi- 
tions about smoothness/flatness of the subsystem boundary. Eq. (3.51) demonstrates that 
for rriL — > 0, AS'free takes the form, 

ASfree = ^ log + 7free (3.54) 

The constant term 7free receives a contribution from modes with energy w ~ 1/L and hence 
cannot be extracted from the present calculation. In fact, given our regularization of the 
trace (3.47), Eq. (3.50) is precisely the difference between the Renyi entropies of a theory 
with mass m <^ 1/L and a theory with mass A <^ l/i^, so the constant 7free cancels out 
in the result. Similarly, with the regularization (3.32), our result (3.42) is the difference 
between the Renyi entropy of a system with spontaneously broken continuous symmetry in 



^ In dimension d — 3, the area law coefficient receives an additional singular contribution, resulting in a 
further correction to the entanglement entropy SS ~ m^^log(ma).^ However, this correction is negligible 
in the regime mL <^ 1, which is of main interest to the present work. 
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a zero external field and the Renyi entropy of a system in a field, which induces a mass gap 
A satisfying Atower -C A <^ 1/L. This is precisely the quantity we considered in Eq. (2.35) 
of section II, which fully agrees with our replica method result (3.42). Moreover, since the 
system in a field is just described by an A'^ — 1 component free scalar field, we may combine 
Eqs. (1.5), (3.54), (3.42) to obtain, 

7ord = (iV - 1) (7free " ^ log TT - + log l^^"^! (3.55) 

Thus, the geometric constant 7ord in the ordered state is directly expressed in terms of 7free- 
We have verified the scahng form, Eq. (3.54), for a lattice regularized free massive bosonic 
theory in d — 2 using the correlation matrix method of Ref. 14, and computed the 7fi.ee term 
in 5*2 as a function of the ratio i/L. The results are shown in Fig. 1. We take the mass 
m = 10~^, the total system as a torus of size L x L and the subsystem as an £ x L cylinder. 
We study systems with L up to 100 sites to extract 7free- 

Our result (3.54) may appear to contradict existing calculations in a free bosonic the- 
ory (see e.g. the review Ref. 26), which find the entanglement entropy to remain fi- 
nite in the massless limit in dimension d > 1, while in = 1 a very weak divergence 
A^* ~ I log log 1/ (m£) is observed. However, the geometry studied in Refs. 26,27 involves 
a subsystem of finite size i embedded in a system whose size L is taken to infinity from the 
outset. The limits m — > and L — > oo do not commute; our result (3.54) corresponds to 
taking m — > while holding the ratio ^/L-fixed. 

IV. CONCLUSION 

In this paper we have demonstrated the presence of logarithmic corrections to the en- 
tanglement entropy in systems with spontaneous breaking of continuous symmetry. Such 
corrections have been recently observed in Monte-Carlo simulations.^'^ Our result, Eq. (1.5) 
gives the coefficient of the logarithmic divergence to be 6 = 1 in the case of a Heisenberg 
antiferromagnet (A^ = 3) in two spatial dimensions. Presently, the Monte-Carlo simulations 
give b — 0.74±0.02.^^ We beheve that the difference between our exact result and the Monte- 
Carlo comes from the difficulty in extracting a subleading correction to the Renyi entropy 
in systems of relatively small size (up to 20 x 20 lattice sites) studied in Ref. 13. In addition 
to the coefficient of the logarithmic correction, we also give an expression, Eq. (3.55), for 
the finite geometric constant 7ord in Eq- (1-5) in terms of a related constant 7free in a free 
massive bosonic theory, see Eq. (3.54). We evaluate the constant 7free as a function of the 
subsystem to total system size ratio i/L using the correlation matrix method of Ref. 14 
(Fig. 1). Wc would like to point out that the behavior of the geometric constant 7ord has 
very recently been investigated by Monte-Carlo simulations of the Heisenberg modeP^ - the 
shape of 7free(^/-^) (and hence of 7ord(^/-^)) in Fig. 1 appears to be at least in qualitative 
agreement with the Monte-Carlo results of Ref. 25. 

A curious byproduct of our work is that the Renyi entropy of a free bosonic theory diverges 
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as the boson mass m — > 0, see Eq. (3.54). This is a consequence of the shift symmetry of 
the free masslcss theory 0^0 + const, (here we assume that the boundary conditions on 
the field do not break this symmetry). If the field (j) is non-compact, the symmetry group is 
M. When one computes an expectation value of an operator O which is invariant under the 
above symmetry, (O) = / D(f)Oe~^^'^y J D(j)e~^^^\ the infinite factor of the group volume 
cancels between the numerator and denominator. However, when one computes the Renyi 
entropy, Sn — ^og{Zn/Z"'), the group volume no longer cancels - this is, essentially, the 
origin of the m — > divergence (3.54). If the field is compact, then the m — > divergence 
is cut-off by the tower of states energy spacing, giving our main result (1.5). Note that the 
divergence of the Renyi entropy in the free non-compact scalar theory (3.54) is present in 
any dimension, including d — 1. Thus, the famous result for the entanglement entropy of 
a free boson in one dimension S = |logL/a implicitly assumes compactness of the boson 
field. The compactification radius is secretly hidden in the short-distance cut-off a. The 
m — )■ divergence of Eq. (3.54) naively appears to contradict previous calculations^^, which 
have found the entanglement entropy of a free non-compact bosonic theory to remain finite 
in the massless limit if the dimension d > 1 (a weak loglog(l/(m£)) divergence was found 
in d— 1^^). However, these calculations have taken the total system size L to infinity first, 
and then m — > while keeping the subsystem size i fixed. This order of limits results in 
the entanglement entropy having the same behavior as in a CFT. On the other hand, in 
Eq. (3.54) wc arc considering a different order of limits, where £/L is kept fixed and m is 
taken to zero. Note that in the case of a state with spontaneous symmetry breaking and no 
external symmetry breaking field we expect that taking L ^ oo with £ fixed simply results 
in a replacement of L by £ in our result (1.5). 

The present paper has mainly focused on the case when the subsystem boundary is smooth 
in d — 2 and flat in d — 3. When this condition is violated, the entanglement entropy is 
known to acquire additional logarithmic contributions. For instance, in the case d — 2, 
when the subsystem boundary has corners, the entanglement entropy of a CFT receives a 
contribution 



Here the sum is over the corners of the boundary and the coefficient bcomifi) depends on 
the corner angle (pi. Such a corner contribution has been found in a free massive bosonic 
theory,® 



where an important condition m » 1/L is assumed. For instance, for a 90° degree angle one 
obtains, fccorn,n=i(7r/2) PS —0.012 for the entanglement entropy proper and &corn,n=2(7r/2) a; 
—0.0062 for the Renyi entropy with n — 2. In a CFT perturbed by a relevant operator 
which produces a mass gap m, the expression (4.2) crosses over to Eq. (4.1) once m <^ 1/L. 
However, as we discussed above, in a free bosonic theory with m — > there is an additional 




(4.1) 




m > 1/L 



(4.2) 
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contribution to the entanglement entropy (3.54), so we expect to logarithmic accuracy^ 

= + E ^-rnlV'O log ^ + ^ log ' ^ « i (4.3) 

i 

Note that the corner coefficients 6corn(</') are controlled by ultra-violet physics and so are 
the same in Eqs. (4.2) and (4.3). The reason for the simple addition of the contributions 
(4.1) and (3.54) to the entanglement entropy is the separation of energy scales from which 
these originate; the contribution (4.1) comes the energy range 1/L ^ w ^ 1/a and the 
contribution (3.54) from m ^ ^ 1/L. As before, for a system with spontaneous breaking 
of continuous symmetry the m — )■ divergence of Eq. (4.3) is cut-off once m becomes of the 
order of tower of states energy spacing Atower ~ —Td so 

Sn = + (AT - 1) J] 6,„rn(^i) log - + log(p,L^- Vc) (4.4) 

Qj a ^ 

In addition to the Renyi entropy, in the present paper we have also studied the entan- 
glement spectrum of systems with spontaneous breaking of continuous symmetry. We have 
shown that the "low-energy" part of the entanglement spectrum takes the same "tower of 
states" form as the physical spectrum of the system. The level spacing of this universal part 
of the entanglement spectrum scales with the system size as c/(psL'^-'^). We hope that this 
finding will be verified by exact diagonalization on small systems or DMRG calculations on 
multi-leg ladders. 
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^ As above, for Eq. (4.3) to be valid, it is important that the to — > Hmit is taken with the ratio £/L - fixed; 
if the total system size L is taken to infinity first and then m — at fixed £ then the CFT expression 
(4.1) is, indeed, obtained as found in Ref. 8. 
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Appendix A: Explicit calculation of Renyi entropy for N = 2 at T = 



In this appendix we show that Eq. (3.43), indeed, remains correct at T = by an exphcit 
calculation in the case N = 2. Here, we may use the angular representation 

fi — (cos 0, sin 0) (Al) 

with the action, 

j d''xdT{d^<Pf (A2) 

The compactness of the order parameter manifold is now expressed through the identification 
(f) 27r. 

Since the theory (A2) is free, we can simply repeat our previous calculation of the entan- 
glement entropy. The only additional complication is that we should restore the periodicity 
of the variable 0. This is accomplished by summing over all winding number configura- 
tions of 0. More precisely, each field configuration on the n-sheeted Riemann surface is 
characterized by the winding numbers r^, with 

/■(fe+l)/3 

/ dr dr(p{T, x) = 27rrfc, 0<x<i (A3) 

and 

/ dr dr(f>{T, x) = 27r V Tfe, £ < x < L (A4) 

^0 k 

Here we are assuming that the field 0(r, x) has no vortices in the (r, x) plane, also known as 
phase-slips. Such phase-slips have an action proportional to L'^^^ and arc, hence, exponen- 
tially suppressed.^ Ordinary spatial vortices are also exponentially suppressed as T ^ 0. 
Any field in the sector with a given set of winding numbers may be written as 

0(x) = 0r(x) + 54){x) (A5) 

Here, 0^ is a reference field carrying winding numbers {vk} and 5(t){x) is a field with all 
winding numbers equal to zero. We choose 0^ to satisfy 

aVr- = (A6) 
We will discuss the solution to this equation shortly. Then 

S[4>\ = S[4>r\ + S[64>\ (A7) 



The phase-slips are also often suppressed by lattice symmetries. 
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Thus, we may perform the integral over 5(f), to obtain 



= Zir'^Wn (A8) 

with 

W^ = Y^e-^^^^^ (A9) 

r 

We have already computed the partition function in the zero winding number sector - 
it is given by Eq. (3.41). It remains to compute S'[0r]- 

Let us solve Eq. (A6) in a sector with winding numbers {vk}- We take to be 

independent of the variables along the subsystem boundary and focus on the (x, r) plain. 
We again introduce fields 0^ via Eq. (3.9). As before the fields are glued along the branch- 
cuts aX. T = k(3 using Eq. (3.10). However, to introduce the winding numbers instead of 
Eq. (3.11) we now impose, 

0fe(r) = 0fc+i(r - ^) + 27rrfe (AlO) 

As before, we work in the limit /3 3> such that the r variable in each sheet runs from 
(— oo, oo). The solution to the Laplace equation with the above boundary conditions is then 
given by Eq. (3.17). In the present case, we may set C = (a finite C can be absorbed into 
50). Using the expansion (3.19) and imposing boundary conditions (3.10), (AlO) we find, 

Cfe = -4-1 (All) 

and 

4-— + (A12) 

Note that due to the relation (All), the constraint (3.18) is automatically satisfied. 

We are now ready to calculate the action S[(f)r] in each winding number sector. We may 
write. 



S[<t>r\ = / dr / d^{d,<j>f = d^{d,<t>,f (A13) 



Here we've broken up the integral into contributions from each sheet. We now integrate by 
parts and use Eq. (A6) to reduce Eq. (A13) to an integral over the boundary of each sheet. 



P L r 
^ k -'^ 

^ dx((0fe9,0fc)(x,^/2) - (0fe9,(/)fe)(x, -^/2)) (A14) 

- / ((0fc5,0fc)(x,O+) - (0fc5,0fc)(x,O-)) 1 (A15) 



k 

2 

k 
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The contribution (A15) comes from the cuts at r = /c/3 and vanishes by Eq. (3.10). On the 
other hand, using Eq. (AlO), the contribution (A14) becomes 



S[<l>r] = ^^^(27rr,) J dxdrM^,P/2) ^ ^^J^^^nn) (-27rc+) « ^ 5^(27rr,) 

k ^ k k 

(A16) 

where we've used Eqs. (3.19),(A12) and dropped a correction of order L//3 in the last step. 
Therefore, the winding number factor Wn in Eq. (A9) satisfies, 

Wn « i/"(e-2"'''^^*^) (A17) 

with = q^'^ - the Jacobi theta function. Observe that 



^ = l + 0(L//3) (A18) 

Hence, the winding number contribution does not to leading order modify our previous result 
for the entanglement entropy Eq. (3.43) for T <^ 1/L. On the other hand, the partition 
function on the n-sheeted Riemann surface is modified by the inclusion of Wn, 

log Zn ""^^ ^ {n/3A + log(np,LV/3) - (n - 1) log(2/3A) + 2nlogi/(exp(-27rV.^''r)) + log(27r)) 

(A19) 

In the limit, T ^ Atowcr ~ ^7?' ~^ ^ partition function (A19) reduces 

to our previous result (3.41) as expected. However, using the identity i/(exp(— a)) = 
■^/tt/q; i/(exp(— tt^/q;)), we may also rewrite Eq. (A19) as, 

logZ„ = i (n/3A - (n - 1) log(47rp,L'^A) + logn + 2nlog z/(exp(-l/(2p,L''T)))) (A20) 

so that for T < Atower, 

log Zn-^^ {npA - (n - 1) log(47rp,L'^A) + log n) (A21) 

Thus, we see that the singular logarithmic dependence of (3.41) in the limit T ^ disappears 
once the sum over winding numbers is performed. 
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